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Abstract. The widely accepted assumption is that metals detected in the spectra of a few cool helium-rich white dwarfs cannot 
K^ . be of primordial origin and therefore must be accreted from the interstellar medium. However, the observed abundances of 

Qs ' hydrogen are much too low to be compatible with the high accretion rates inferred from metal accretion. Hydrogen accretion is 

r\J I therefore suppressed compared to metal accretion. The hypothesis most widely discussed as cause for this "hydrogen screening" 

^SJ ■ is the propeller mechanism: Metals are accreted in the form of grains onto a slowly rotating, weakly magnetized white dwarf, 

^^ ' whereas ionized hydrogen is repelled at the Alfven radius. We have obtained circular polarization spectra of the helium-rich 

O ! white dwarfs GD 40 (WD0300-013) and L745-46A (WD0738-172) — which both show strong metal lines as well as hydrogen 

\l ■ — in order to search for signatures of a weak magnetic field. The magnetic field strengths necessary for the propeller mechanism 

to work in these stars are at least 144000 G and 3000 G, respectively. Whereas L745-46A might have a magnetic field of about 
-6900 G no magnetic field could be found with an upper limit for the field strength of 4000 G (with 99% confidence) for 
Cii' GD40. 
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1. Introduction white dwarf may experience highly variable accretion rates 

during its life. Most of the time (about 5-10^ years assumed 

Some of the helium-rich white dwarfs at the cool end of the -^ ^j^^ ^^.^„gjy idealized scenario) the white dwarf will accrete 

H white dwarf sequence show evidence for metal lines in their ^^ ^^^^ j^^ ^^^^^ (5 . ^q-20 ^g/yr) which lead to metal abun- 

spectra. Since radiation forces are not strong enough to com- ^^^^^^ -^ ^j^^ photosphere well below the detection limit. Only 

pete with gravity for temperatures below 40000 K it is expected ^^^-^^ ^^^ ^j^^^^j^ ^^^^^^ encounters with denser patches in the 

that metals sink down on time scales which are short compared interstellar medium (typical density 10 particles cm^ and tem- 

to the cooling age. If nevertheless metals are observed in the pe^ature 100 K), metals are accreted with rates high enough to 

atmospheres, as in the DZ and related spectral types, they must le^^ to detectable features in the spectra, 
have come from outside the star. The most popular mechanism 
for this is accretion from interstellar matter. 

Based on earlier studies (Koester [79761 Wesemael [M3 Whereas this scenario can account for the observed metal 

Vauclair et al. 1979) concer ning accretion proces ses onto cool abundances in helium-rich white dwarfs, the observed abun- 

white dwarfs, Dupuis et al. JT^ [1993^1 rT993b !) developed a dances of hydrogen are much too low to be compatible with 

theoretical model which predicts element abundances in cool the acc retion rates inferred from metal accretion (Dupuis et al. 

white dwarf atmospheres under the assumption that elements T^^^'- I" ^ recent study Wolff et al. ("MS i analyzed the UV 

are accreted in solar element proportions from the interstellar sP^ctra of ten DZ and DBZ stars which form a homogenous 

medium. They followed the suggestion of Wesemael (1979) data set and concluded that the observed hydrogen abundance 

that due to the patchy structure of the interstellar medium a is more than two orders of magnitude to low. Furthermore DZ 

white dwarfs are not found in or close to known dense clouds 



Send ojfprint requests to: S. Friedrich (Aannestad et al. I1993> . One possible explanation is that hy- 

* Based on observations collected at the European Southern drogen accretion is reduced compared to metal accretion. The 

Observatory, Paranal, Chile (ESQ Programme 66.D-0541) mechanism most widely discussed as reason for this "hydrogen 
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screening" is the propeller mechanism of Illarionov & Sunyaev 
(fT975t. 



2. The propeller mechanism 

The propeller mechanism was originally developed for ac- 
cretion onto neutron stars (Illarionov & Sunyaev 1975) and 
adopted to white dwarfs by Wesemael & Truran (1982 1. 
Wesemael & Truran proposed that metals are accreted in the 
form of grains onto a slowly rotating, weakly magnetized white 
dwarf, whereas ionized hydrogen is repelled at the Alfven ra- 
dius. Fig. 1 in Wesemael & Truran shows the domain of al- 
lowed Alfven radii as a function of luminosity (or temperature) 
of the central star. The domain is restricted to higher tempera- 
tures by setting the dust evaporation time scale at the Alfven 
radius equal to the free-fall time from the accretion radius, 
to lower temperatures by equating the ionization radius to the 
Alfven radius. For typical white dwarfs (M = O.6M0) and ac- 
cretion parameters (M = lO^'^Mo/yr, v - 50km/s) this region 
is roughly between T^ff = 7500 K and T^ff = 15000 K. The 
lower boundary can be bypassed if accretion from a warm ISM 
phase with temperature T ^ 10000 K and density n sj 1 cm"^ is 
assumed. In such a phase, the hydrogen would be ionized, but 
the metals may still be partially locked up in grains (Alcock 
& Illarionov .1980) . Given the effective temperature of a star 
the minimum Alfven radius, which in turn corresponds to a 
minimum magnetic field strength necessary to let the propeller 
mechanism work, can be calculated. 

To our knowledge it has never been tested whether in 
helium-rich metal-line white dwarfs the necessary conditions 
for rotation, magnetic field, and strength of hydrogen-ionizing 
radiation are indeed fulfilled to make this mechanism operate 
over the entire range of Tg^ from 7500 K to 15000 K. 

3. Observation and data reduction 

With effective temperatures of 15150 K and 7500 K, re- 
spectively, GD40 (WD0300-013, V=15.56) and L745-46A 
(WD0738-172, V=13.04) mark the upper and lower boundary 
of the temperature regime of the propeller mechanism. Both 
stars are well studied and atmospheric parameters as well as 
element abundances are known with sufficient accuracy. Based 
on the accretion-diffusion model of Dupuis et al. Cl993bj the 
hydrogen accretion rates necessary to provide the observed hy- 
drogen abundances are about three or four orders of magnitude 
lower than the average accretion rate which is assumed to ac- 



count for the observed metal abundances (Friedrich et al. ll999l 
Koester & WolfF2000"l. Therefore, hydrogen accretion must be 
efficiently reduced even onto L745-46A which shows one of 
the highest hydrogen abundances observed in cool heUum-rich 
white dwarfs. 

For the above reasons both stars are ideal candidates to 
search for magnetic fields driving the propeller mechanism. 
The VLT with FORSl and PMOS offers for the first time 
the possibility to obtain circular polarization spectra with high 
enough signal-to-noise ratios at acceptable exposure times to 
search for these magnetic fields as it has been demonstrated by 
Aznar Cuadrado et al. (2004) for white dwarfs. 



Flux and circular polarization spectra were obtained with 
the VLT UTl and PMOS in December 2000 in service mode. 
Grism 600R (dispersion lA/pix) was used with order sorting 
filter GG435 resulting in a spectral range of about 5300-7350A. 
Several individual exposures were taken for each star. After 
each exposure the retarder plate was rotated by 90 degrees in 
order to account for improper flat field corrections, instrumen- 
tal polarizations and depolarization. The resulting total expo- 
sure times for GD40 and L745-46A were 3.3 hours and 2.7 
hours. The signal-to-noise ratio in the vicinity of Ha and the 
helium lines of GD40 amounts to 190 and 260, respectively, 
and to 280 for Ha of L745-46A. Data reduction was done with 
standard IRAF routines. Circular polarization was first calcu- 
lated for each individual retarder position from the difference 
of the double spectrum, produced by the Wollaston prism, di- 
vided by the sum of the double spectrum. Polarization spectra, 
corrected for instrumental influences, were then derived from 
the difference of two consecutive circular polarization spectra 
with different quarter wave plate positions. Subsequently these 
polarization spectra were averaged to get a mean circular po- 
larization spectrum. The flux spectrum is given by the sum of 
both spectra. All individual flux spectra were averaged to get a 
mean flux spectrum. 

In addition a star in the field of view of L745-46A, which 
we assume to be unpolarized, was observed. It served as a sup- 
plementary check for instrumental introduced polarization. The 
mean circular polarization in a spectral range between 5500A 
and 6800A amounts to 0.01%, 0.02%, and 0.03% for the field 
star, GD 40, and L745-46A, respectively, one order of mag- 
nitude lower than the expected values for kilogauss magnetic 
fields. 

4. Determination of magnetic fields 

4.1. Weak-field approximation 

According to the theory of line formation in a weak magnetic 
field the splitting of a spectral line is proportional to the mean 
longitudinal field B\, i.e. the component of the magnetic field 
along the line of site averaged over the visible stellar surface. 
Provided that the Zeeman splitting of a spectral line is small 
compared to intrinsic (thermal and pressure) broadening (e.g. 
Angel & Landstreet 1970) the amount of polarization can be 
determined by the weak-field approximation, which is valid 
even in the presence of instrumental broadening, but is not gen- 
erally correct in the case of rotationally broadening (Landstreet 
[19821 : 
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with elAimieC^ ^ 4.6710"''' A"' G"', and geff the effective 
Lande factor which equals 1 for hydrogen Balmer lines (e.g. 
Bagnulo et al. 2002'! and citations therein). For the Hei line at 
5875A geff = 1-16 (Leone et al.'2000'). 

In our case the weak-field approximation is justified be- 
cause intrinsic broadening is much larger than the separation of 
line components due to Zeeman splitting which is of the order 
of 0. 1 A for field strengths of 10 kG. For higher field strengths 
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which are called for by the propeller mechanism for GD40 
the intrinsic broadening is smaller than the magnetic broaden- 
ing and the weak-field approximation reduces to a order-of- 
magnitude estimation (Landstreet .198^ . 

Due to the lower signal-to-noise ratio at the Ha line of 
GD40 the synthetic flux spectrum I{A) from Section 5.1.1. 
was used to calculate V/I. For L745-46A the S/N-ratio was 
sufficiently high to use the observed flux spectrum. dl/dA is 
approximated by the average of (/,+i - /,)/(/l,+i - Aj) and 
(/,■ - /,_i)/(/l,- - /i,_i) for a given wavelength Ai. 

4.2. Fitting of the magnetic field 

In order to determine the mean longitudinal component of the 
magnetic field the observed circular polarization was compared 
to the predictions of the weak-field approximation in an inter- 
val of +20A around Ho- and the He i lines at 5875A and 6678A. 
The best-fit for B\, the only free parameter, was found by ax^ 
minimization procedure. If no magnetic field were present, all 
deviations from zero would be due to noise. This can be ex- 
pressed by the standard deviation over the respective intervals 
around the respective lines. If the best-fit magnetic field is close 
to zero the reduced x^ should be close to unity. However, if a 
magnetic field is present some deviations from zero are not due 
to noise and the reduced x^ is smaller than 1 . Following Press 
et al. ([^86 1 the statistical error is determined from the rms 
deviation of the observed circular polarization from the best- 
fit model. The Icr (68.3%) confidence range for a degree of 
freedom of 1 is the interval of B\ where the deviation from the 
minimum is Ax^ = 1 (6.63 for the 99% confidence range). This 
is a purely statistical error and does not account for systematic 
errors. For details see Aznar Cuadrado et al. C2004J. 
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Fig. 1. Observed flux spectrum (dashed) of GD40 at the po- 
sition of Ha with a non-magnetic model spectrum (solid, 
^eff ~ 15150, logg - 8) over-plotted. The hydrogen abun- 
dance needed to reproduce the line profile of Ha amounts to 
5 ■ 10"^ relative to helium 



determined an equivalent width of 1 A, we measured 0.2 A only. 
We attribute the difference in the equivalent widths to the much 
better quality of our spectra compared to those of Greenstein 
& Liebert. As a consequence the hydrogen abundance turned 
out to be a factor of 10 below the value previously determined 
(Friedrich et al.'1999). In order to estimate the hydrogen abun- 
dance, we took the best fit, non-magnetic model atmosphere 
with Tgff = 15150, logg = 8 from Friedrich et al. iT999t 
and varied the hydrogen abundance while keeping all other el- 
ement abundances and the atmospheric parameters fixed. This 
resulted in a new hydrogen abundance of 5 ■ 10"^ relative to 
helium for GD 40. 



4.3. Magnetic field strengths for the propeller 
mechanism 

Minimum magnetic field strengths necessary for the propeller 
mechanism to work were calculated according to Wesemael 
& Truran <1982> . The logarithm of the Alven radius was es- 
timated from their Fig. 1. Assuming a white dwarf mass of 
O.6M0, a radius of O.O13/?0, an accretion rate of 10"'''Mo/yr 
the corresponding minimum magnetic field strengths from their 
Eq. 4a are 144000+^^;j[J° G and 3000+500 G for GD40 and 
L745-46A, respectively. The error accounts for a reading er- 
ror of 1mm in their Fig. 1 . Further errors are introduced by the 
stellar radius which enters the formula with third power and the 
accretion rate. Lower radii and higher accretion rates would in- 
crease the minimum magnetic field strength. 



5. Results 

5.1. GD40 

5.1.1. Hydrogen abundance 

The flux spectrum of GD 40 clearly shows an Ha absorption 
line which confirms the detection of this line by Greenstein & 
Liebert (1990 Fig.^. In disagreement to these authors, who 



5.1.2. Circular polarization 

The inspection of the individual polarization spectra of GD 40 
did not show any obvious features or variations in the polariza- 
tion which cannot be attributed to noise. The magnetic field 
strengths calculated with the theoretical weak-field approxi- 
mation are below the Icr error In addition different magnetic 
field orientations are found (indicated by positive and nega- 
tive signs of Bi) for Ha and the two He: lines. For the Ha 
line the field strength and the respective 68.3%/99% confi- 
dence ranges amount to (2747 ± 6892/15000) G, and for the 
two He I lines at 5875A and 6678A to (-1683 + 1528/3936) G 
and (-716 + 1848/4758) G, respectively. A fit to all three 
lines results in a magnetic field strength of 1 131 G and a 99% 
confidence level of + 2973 G. For comparison the derived 
field strengfli for the field star amounts to -438 + 318/819 G 
(68.3%/99% confidence range). Thus we conclude that GD40 
does not posses a magnetic field with an upper limit of about 
4000 G. 

Fig. 12 shows the observed circular polarization spectra of 
GD 40 at the Ha line and the two He i lines together with the 
respective polarization spectra with the magnetic field strengths 
predicted by the weak-field approximation and the maximum 
field strengths within the 68.3% and 99% confidence ranges. 
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5.2. L745-46A 

Fig-E]shows the circular polarization at Ha, which is the only 
visible line in the observed spectral range, together with the 
polarization spectrum for the best-fit magnetic field strength of 
-6900 G. The 68.3% and 99% confidence ranges amount to 
+2100 G and ±5550 G. Again the polarization spectra with the 
maximum value from this field strength and the 68.3% and 99% 
confidence range are also plotted. Although the field strength 
clearly exceeds the statistical Icr error one has to be cautious, 
because the result solely depends on the vicinity of the Ha line. 
However, if this field could be confirmed, it would be strong 
enough to prevent hydrogen from accreting according to the 
prediction from the propeller mechanism. 

6. Conclusion 

We have obtained high signal-to-noise circular polarization 
spectra of the helium-rich white dwarfs GD 40 and LP745-46A. 
With errors as low as +0.5% and ±0.2% these spectra allow for 
the first time to search for the weak magnetic fields necessary 
to drive the so called propeller mechanism which is assumed to 
be responsible for the low observed hydrogen-to-metal ratios 
found in some cool helium-rich white dwarfs. 

For L745-46A the predicted magnetic field strength from 
the weak field approximation amounts to -6900 G with a 
formal Icr error of 2100 G and a 99% confidence range of 
±5550 G. This magnetic field strength exceeds the minimum 
field strength of 3000 G necessary for the propeller mechanism 
to work. However, one should keep in mind, that this result is 
based on the Ha line only. If the detection is confirmed, this 
means for the first time an indication that magnetic fields may 
play a role in the accretion of hydrogen and metals onto cool 
helium-rich white dwarfs. 

However, we could not detect signatures of a magnetic field 
in the circular polarization or flux spectrum of GD40. This 
does not necessarily mean, that no magnetic field is present, 
because if we are looking on the magnetic equator of a mag- 
netic dipole (/ - 90°), the components of the magnetic field 
along the line of sight completely cancel and no circular po- 
larization can be detected. According to Wesemael & Truran 
(1982) a field strength of about 144000 G is needed for GD 40 
to let the propeller mechanism work. If we now assume that 
we could detect a field strength of 4000 G (derived upper limit 
field strength from the fit to all three lines with a confidence 
of 99%) we could estimate an inclination angle of greater than 
84° for which the longitudinal field strength becomes to low to 
be detected (Brown et al. 1977 1. Therefore the chance to miss 
a magnetic field of 144000 G is only about 6%. 

From our observations we must therefore conclude that at 
least for GD 40 other mechanisms than magnetic fields prevent 
hydrogen from accreting. Sion et al. ( 1990 l argue that hydrogen 
deficient stars like GD 40, which do not lie close to any known 
interstellar cloud (Aannestad & Sion "1985'l, have accreted, 
but not from interstellar clouds. They proposed that white 
dwarfs might experience accretion in different environments 
with some white dwarfs being accreting material from interstel- 
lar clouds and some from volatile-free asteroidal or dust mate- 



rial or chemically differentiated planetary material. Accretion 
of comets from an Oort-like cloud as proposed by Alcock et 
al. (1986 1, however, failed to explain the observed Ca/H ra- 
tios and is ruled out by Zuckerman et al. t2003j as a source of 
metals at least for DAZ white dwarfs. In view of the sugges- 
tion that the DAZ phenomenon might be associated with the 
existence of a nearby unseen companion (Zuckerman & Reid 
1998), Zuckerman et al. (2003 1 also speculate about Jupiter- 
size substellar companions of DA and DB white dwarfs, that 
might generate detectable metal-line absorption. One problem 
with this model is lifting material off the secondary. Another 
problem might be the lack of observed eclipses, which would 
be expected, if a substantial number of all DAs and DBs are 
orbited by such substellar companions with semi-major axis of 
a few solar radii, which are typical for red dwarf/white dwarf 
pairs. A dusty circumstellar disk might be also the reason for 
a hydrogen depleted environment as mentioned in Aannestad 
et al. ( 1I993> . The maintenance of such a disk over the long 
lifetime of a cool white dwarf, however, would be difficult to 
understand. 

From the determination of hydrogen and metal abundances 
in helium-rich white dwarfs one must conclude that predomi- 
nantly metals are accreted. However, the reasons for this selec- 
tive accretion process are not yet really understood. It seems 
that also other mechanism than magnetic fields must be con- 
sidered. 
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Fig. 2. Observed (gray) circular polarization spectrum of 
GD 40 together with the respective polarization spectra (bold) 
predicted by the weak-field approximation for field strengths of 
2747 G, 9639 G, and 17747 G for Ha (top), - 1683 G, -321 1 G, 
and -5619 G for the He i fine at 5875A (middle), and -716 G, 
-2564 G, and -5474 G for the He: line at 6678A (bottom). 
For each spectral line the first field strength denotes the best 
fit value, the second and third value the best fit value plus the 
68.3% and 99% confidence range, respectively 
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Fig. 3. Observed (gray) circular polarization spectrum of L745- 
46A at Ha together with the polarization spectra (bold) deter- 
mined by the weak-field approximation for field strengths of 
-6900 G, -9000 G, and -12450 G. The first field strength de- 
notes the best fit value, the second and third value the the best fit 
value plus the 68.3% and 99% confidence range, respectively 



